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Abstract 

We record what is known about the closed forms for various Bessel function moments 
arising in quantum field theory, condensed matter theory and other parts of mathematical 
physics. More generally, we develop formulae for integrals of products of six or fewer Bessel 
functions. In consequence, we are able to discover and prove closed forms for Cn^k '■= 
/q°° t^Xo (t) dt with integers n = 1,2,3,4 and k > 0, obtaining new results for the even 
moments C3 2k and C4 2k ■ We also derive new closed forms for tlie odd nionients 572,2/^+1 • — 
f^t^''+^lQ(t)K^-\t)dt with n = 3,4 and for tn,2k+i ■= t^''+^ll{t) K^-'^{t)dt with 
n = 5, relating the latter to Green functions on hexagonal, diamond and cubic lattices. 
We conjecture the values of s^^2k+i-, make substantial progress on the evaluation of c^^2k+ii 
SQ,2k+i and i6,2fc+i and report more hmited progress regarding c^^2k, C6,2fc+i and CQ^2k- 

In the process, we obtain 8 conjectural evaluations, each of which has been checked 
to 1200 decimal places. One of these lies deep in 4-dimensional quantum field theory and 
two are probably provable by delicate combinatorics. There remains a hard core of five 
conjectures whose proofs would be most instructive, to mathematicians and physicists 
alike. 
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1 Introduction 



We give representations of the vacuum [35l [56] and sunrise [30l HOj E] diagrams 




where all the arguments of Vn and all but the last argument of Sn+i are real and positive. 
Here and below Iq, Jq and Kq are the conventional Bessel functions of order zero as in [H 
Chapter 15]. These integrals occur in quantum field theories in two spacetime dimensions, 
where we do not need to regularize ultraviolet divergences. Numbers generated by them 
are expected to occur in the finite parts of integrals from Feynman diagrams in four 
spacetime dimensions. To be concrete, we illustrate V3 and S'4 as follows: 





By casting Bessel's differential equation in the form 

- + -^] -^Ko{at) =t^Ko{at) 
a day da 

and applying the corresponding differential operator to ([1]), we may increase the exponent 
of t in the integrand by steps of 2. But to obtain even powers of t, we need to start with 

/■°° / " \ 
Vniai, ...,an):= (J] ^0(0^*) I dt (3) 

which plays no obvious role in quantum field theory. To evaluate the latter form, we found 
it useful to regard the Fourier transform 

Sn+i{ai, . . . ,an,w) := - I \ f] K q {a jt)\ cos (wt) dt (4) 

V=i J 

as an analogue of ([2]). 

We shall be especially interested in the moments 

c„,fc:= / t%"(t)dt (5) 
Jo 

for integers n > 1 and A; > 0, as studied in [HI |TD] and [261 160] . 
In [To] these moments arose in the study of Ising-type integrals 

_ 1 dxi dx2--- dx 



Jo Jo (cosh xi + • • • + cosh Xn)^~^^ 
2 



which are hnked by 

Cn.k = rrr ^n.k ■ (6) 

In \26\ [60] it is proven that for fixed n these moments satisfy a hnear recursion for 
which a simple algorithm exists with coefficients polynomial in k. For example, for n = 1 
and 2 one easily obtains the closed forms 

\ 2 J 4r(| + i) ^ ^ 

and for n = 3 and 4 we obtain the recursions 

(A: + - 2(5A;2 + 20A; + 21)c3,fc+2 + 9c3,fe+4 = (8) 

(A: + l)5c4,fc-4(A: + 2)(5A:2 + 20A: + 23)c4,fc+2 + 64(A: + 3)c4,fc+4 = 0. (9) 

These recursion formulae may be written quite compactly as 

M 

J2i-'^yPn,i{k + i+^)Cn,k+2i=0 (10) 
i=0 

where M = [{n + l)/2j . For instance, for n = 5 and 6, we have 



P5 = 353;^ + 42x2 + 3 v%^\[x) = x(56x'' + ll2x^ + 24) 
P5^2(x) = 259x2 + 104 P6,2{^) = x(784x2 + 944) 
P5,3{x) = 225 pe,3{x) = 2304x. 











(11) 



The same recursions apply to the moments 

/•oo 

Sn,k := I tHo{t)K^-\t)dt (12) 
for integers n > 3 and A; > and to 

poo 

tn,k ■■= I t''liit)K^-\t)dt (13) 

for integers n > 5 and A; > 0. 



2 Two Bessel functions 

The transforms 

S2{a,w) := / tKo{at)Jo{wt)dt = ^ ^ (14) 

Jo + 

82(0, w) := - I Ko{at) cos{wt) dt = — . „ (15) 
give Vi{a) = l/o? and Vi{a)/T: = l/(2a) at = 0. The distributions 

poo 

/ vUvh)Jo{vt2)dv = 25{tl-tl) (16) 
Jo 

2 

- / cos(T;ti) cos(z;t2) dt; = 5{ti + 12) + 5{ti - t2) (17) 
v" Jo 



then lead to the evaluations 

fOO POO 



V2{a,b) := / tKo{at)Ko{bt)dt = / wS2{a,w)S2{b,w) dw 
Jo Jo 
log(a/6) 



/"OO /"CxD 

F2(a,6) := / Koiat)Ko{bt)dt = 2TT / ^2(0, w):S2(6, u;) 
Jo Jo 



(18) 



K ^ (19) 



a + 6 \o + 6 
with a complete elliptic integral of the first kind, 

Kik):= _ dcf>, (20) 

JO V 1 — K"' sm 

appearing in (jl9p and a limit intended in (jlSp when a = b. We shall need to refer to 
the complementary integral K'(A;) := K.{k'), with fc' := \/l — /c^. In the case a > b, this 
provides a compact alternative form of p^ . 



y2(a,6) = ^K'(6/a), (21) 
obtained by the Landen j3j transformation in [T, 17.3.29]. 

3 Three Bessel functions 

We follow Kallen [511 [52] by constructing S3{a,b,w) from its discontinuity across the cut 
in the w'^ plane with branch point at w'^ = — (a + 6)^, obtaining 

Ss{a,b,w)= tKo{at)Ko{bt)Mwt)dt= 7 ' J ^ d^ 

Jo Ja+b V +W 

with a discontinuity 

Ds(a,b,c) = — , 

^y{a + b + c){a -b + c){a + b-c){a-b- c) 

that is completely symmetric in its 3 arguments. The u-integral is easily performed, to 
give 

53(a, b,w)=2 arctanh ^^ ^^ | |" ~ Ds{a, b, iw) (22) 
and in particular the on-shell value [40J 

/•OO 

S3,i = 53(1, l,i) := tIo{t)Kl{t)dt = L_3(l) = — 



3^/3 



where 

X-3("') / 1 1 



, A?,k + lY (3A; + 2) 

n=l A:=0 ' ^ ' 

is the Dirichlet L- function with the real character x-'ik^^ given by the Legendre-Jacobi- 
Kronecker symbol (-D|n) for discriminant D = — 3. 
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Nomenclature. We shall refer to the construction of a Feynman amplitude from its 
discontinuity across a cut as a dispersive calculation. Kramers [53] and Kronig [5^ 
founded this approach in studies of the dispersion of light, in the 1920's. In the 1950's, the 
utility of dispersive methods was recognized in particle physics |45 1 I52 1 [69]. Cutkosky [31] 
turned them into a calculus that became a routine part of the machinery of quantum field 
theory. Barton [20] has given a scholarly and instructive introduction to these techniques. 



3.1 The odd moments 53 2^+1 

By differentiating 

2 arcsin(x/2) 

'53(1,1,12;) = ^ 

a;v4 — 

before setting x = 1, we evaluate 

1^ 1 , d ^ dSsil, l,ix) 



= o ■53,1 
x=l "J 



dx 

and are then able to solve the recursion relation |26l [60] for 53^2^+1 by the closed form 



with integers 

Integer sequence (p^ begins 
1, 3, 15, 93, 639, 4653, 35169, 272835, 2157759, 17319837, 140668065, 1153462995 (25) 
and is recordec0 as entry A2893 of the on-line version of [66] , which gives the recursion 

{k l)^afe+i - (lO/c^ + lO/c + 3)afc "^k^a^-x = (26) 
and the generating function 

Il{2t) = Y,aA-\ . (27) 
fc=o ^ 

We have verified that recursion (j26p reproduces the recursion for (j23p , which has the same 
form as for the odd moments in We note that integers at were encountered in studies 
of cooperative phenomena in crystals [36] and also in studies of matrices [43] with entries 
or 1. In [321 Prop. 2], they are related to enumeration of closed walks in a two-dimensional 
hexagonal lattice. They also appear in an enumeration of Feynman diagrams |6H Table 
2] in quantum chromo dynamics, via the constrained sum [18] 

V- f k\ 
ak= 2^ 

p+q+r=k 



pi q\ rl 



that results from the Taylor expansion of Iq. 

It is notable that Iq provides a generating function for moments of IqKq. In Section 4 
we shall show that Iq generates moments of IqKq. 



^ See http : //www . research . att . com/ ~n j as/ sequences/ A002893 
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3.2 The odd moments C3 2fc+i 

Next, we construct 

roo poo 

Vs{a,b,c) := / t Ko{at)Ko{bt)Ko{ct) dt = / wSs{a,b,w)S2{c,w) dw 
Jo Jo 

_ L3{a, b, c) + ^3(6, c, g) + 13(0, a, b) 
4 

with a dilogarithmic function 



-1)3(0, 6, c) (28) 



L3(a, 6, c) := L12 , ^ . ^ 9Y7T7 — 7 " L12 ' 



(a2 + 62 _ c2)L>3(a, b,c)-lj '\{a^ + b'^- c^)D3{a, b,c) + l 
computed by Davydychev and Tausk [35] . Setting a = 6 = c= l,we obtain 

n 5 



C3,l 



l/3(M.i) = ji-s(2) = iE(i^)"E 

n=0 ^ ^ fc=l 



(6n + kY" 



where the vector of coefficients v = [9, —9, —12, —3, 1] was discovered (and proven) in the 
course of investigation of 3-loop vacuum diagrams in 4 dimensions [29j . 
Similarly, 

C3,3 = L^z{2) - ^ 

may be obtained by suitable differentiations of (j28p . Then higher moments 03^2^+1 with 
k > 1 may be obtained by using ([8]). Because of the mixing of L_3(2) with unity, we were 
unable to write a closed form for their rational coefficients in C3^2fc+i- 



3.3 The even moments C3 



2k 



For even moments, we lack the dispersion relations |20j of quantum field theory and so 
fall back on the general Aufbau 

/oo 
'S'm+i(ai, ■ ■ ■ am,v)Sn+i{bi, ...bn,v + w)dv (29) 
-00 

which follows from the distribution 
1 



vr 



00 



cos(fti) cos((t> + w)t2) dv = 6{ti + t2) cos{wt2) + S{ti — 12) cos{wt2) 



obtained from (fT7|) and the expansion cos[{v+w)t2) = cos(t;t2) cos(i(;t2)— sin(t;t2) s\ii.{wt2) ■ 
In particular, by setting m = n = 1 in (j29p . we evaluate 



vr 



00 



S3{a,b,w) := - I Ko{at)Ko{bt) cos{wt) dt = I S2{a,v)S2{b,v + w) dv 







as the complete elliptic integral 



1 dv K IV 



S3{a,b,w) = - -.= ' ' ^ (30) 

4 V(a2+^;2)(62 + (^ + y,)2) ^^2 + + 5)2 
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and recover the previous result (fT9|) for V2{a, b) = 7rS3(a, b, 0) by setting w = 0. 
Next, by setting m = 2 and n = 1 in ([29]) . we write 



S4{a,b,c,w) := — I Ko{at)KQ{bt)KQ{ct) cos{wt) dt = I S3{a,b,v)S2ic,v + w) dv 



IT 



as an integral over an elliptic integral: 



S,(a, c. = 1 r —^M^^^= d„ (31) 
2 ^(„2 + (a + (,)2)((i- + ,i,)2 + c2) * ' 

and at = obtain 

K 



Vs{a, b, c) := r Ko{at)Ko{bt)K,{ct) dt = vr ^^V^W^ _ ^32) 
Jo Jo V(i^2 + (a + 6)2)(i;2 + c2) 

Remarkably, the integral (j32p may be evaluated by exploiting a more general identity 
given in W.N. Bailey's second paper on infinite integrals involving Bessel functions |17j . 
Without loss of generality, we assume that c > 6 > a > and define 



2c 



2 



Then our result may written as 
2c 



^3(0, b, c) = K(/c_) K{k') + K(A;+) K(A;'_) . (33) 

TT 



We remark that when c > a + b each term in ([33]) is real; otherwise, each is the complex 
conjugate of the other. The form of k± comes from Bailey's conditions k^k^ = a/c and 
k'^k'_ = b/c. 

We proved ([33]) by setting 1/ = p = in Bailey's equation (3.3), to obtain 

1 

I^{at)K^{bt)Ko{ct) dt = -W^{k+)W^{k.) 
with a hyper geometric series 



n=0 



where P is the Legendre function defined in [T, 8.1.2]. Then ([33|) follows from the expan- 
sions 

I^{x) = Io{x) - fiKoix) + 0{fi^) 
W^{k) = 2K{k) - pvrK ^1 - k"^^ + 0{iJ.^) 

where the derivative of Wf^{k) at = is obtained by setting a = b = ^ and z = 1 — k"^ 
in [U 15.3.10]. 
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Specializing (j33p to the case a = b = 1 we obtain 
2c 



vr 



K^{t)Ko{ct) dt = A{2/c) = B{c/2) 



(34) 

(35) 
(36) 



with (I35p coming directly from (1330 . The identity ^(x) = -B(l/j;) may be proven by 
showing that ^(1 — y) and 5(1/(1 — y)) satisfy the same third-order differential equation 
and have Taylor series about y = that agree in their first 3 terms. 
Alternative one may use transformations of the Meijer G-functions 



with the choice of a sum of squares or a product in the functions 
Aix) := K^(^^-^)+K 

/ A /I -A- X — A /I — T. 

B{x) : 



23 




1 


1 


1 


33 




1 


1 


1 




yX^ 


2 


2 


2 


23 




1 


1 


1 


33 




1 


1 


1 






2 


2 


2 



i32 
'33 





1 


1 






2 


2 













^23 
'33 





1 


1 






2 


2 













(37) 
(38) 



in the notation of [63, Vol. 3]. This provides an inversion formula, 8.2.1.14, used in ([3 
and a multiplication formula, 8.2.1.15, used in (1380 . 

Remark. Our forms, ([35]) and ([36]) . for the Bessel integral ([M|) . were tabulated, without 
proof, in identities 2.16.46.4 and 2.16.46.5 of [63] Vol. 2]. Our proof of the more general 
identity (j33p came from following a reference to Bailey's work, given in Section 7.14.2.43 
of the Bateman project [37l Vol. 2]. 



(39) 



Setting c = 1 in (j34j) . we obtain 



/ 3r6(i) 

C3,o = 2 ^3i^3 = 



with the product of K-^ = K(A;3) and K'^ = V^K^ obtained from (j36p . at x = ^, where 
the third singular value |24j 

v3 — 1 / , N 

' = = sm(^/12) 

results. Moreover, Bessel's differential equation yields 



TT /I d \ d 5(c/2) 
C3,2 = - I - + 



2 \ c dc / dc 
upon use of the evaluations 

E(sin(7r/12)) 
E(cos(7r/12)) 



c=l 



r6 (1) 47r522/3 
967r22/3 9r6(i) 



(40) 



6 



+ 



vr 



4i^3 
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of complete elliptic integrals of the second kind, recorded in |24l p. 28] and first found by 
Legendre. Prior to finding this proof, we discovered ()40p in the more palatable form 

1 -1 

C3,2 = g C3,0 - ^ Cg (41) 

by using PSLQ O [121 112] in a manner suggested by previous discoveries in quantum field 
theory, as described in Section 5. 



3.4 Continued fraction 

We recall that 



92«= (42) 

C3,o „ , 9-3^ 



d{l) 



9 • (2iV - 1) 



d{N)-... 

where d{N) := 40A^2 _)_ 2 was derived in p^. Hence, dividing (j4ip by cs^O) we obtain a 
neat continued fraction for 

l-2V39('£j|iy = 9^. (43) 

3.5 Double integrals 

Inspection of [lOj also reveals that we now have evaluated the integrals 

The first integral occurs in a formula for 4c3,o/(3vr) in ^Q\. To evaluate the second integral, 
we note that it is twice the value obtained from the simplex x > y > 0. The transformation 
X = y + z then proves that (f45]) is twice (f44l) . 

In Section 5, we give evaluations of double integrals arising in quantum field theory. 



3.6 Hypergeometric series 



We may also obtain a simple hypergeometric series for cs^o from the Clausen product 
formula |24l p. 178] in the form 



-2K2(sin(a/2))=3F2( 2'^2|2 



1 1 1 



sin^ a 



(46) 



which is valid for 7r/2 > q > 0. Setting a = tt/6, so that sin(a) = 1/2, we recast ([39]) as 

f2n\3 



3.0 g 8 

n=0 



111 
2' 2' 2 
1,1 



(47) 
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Moreover, we conjectured the compact formula 



^'^ 288 ^ 28" (n + 1)2 288 ^ H 2,2 

n=0 V I ; \ ' 



j) 

as an alternative to (j4ip . A proof was provided by Maple, which evaluates this sum in 
terms of and E^. Our evaluation of the latter then shows that ()48p follows from (I40p . 

3.7 Integrals of elliptic integrals 

Setting a = 6= l,i(; = 2 tan and c = 2 sin a in (j32p , we obtain from (j34p the evaluation 

K(sin^) 5(sina) 

/ = du = 

^0 cos e V tan2 6 + sin^ a 2 sm a 

and hence, by trigonometric simplification, we prove the identity 

/ . , \, = K(sin(a/2)) K(cos(a/2)) (49) 

JO V 1 — cos^ a cos^ 

which we had conjectured empirically from interpolating known results, using Maple^s 
MinimalPolynomial. At a = vr/2, identity (j49p reduces to the evaluation 







\V2j 167r 
given in [2U p. 188]. At a = 7r/6, we obtain 



2r-^m=aS^K,Ki. (50) 

V 1 + 3 sin^ 9 



We did not find identity ()49p in the literature. However, John Zucker remarked to us 
that the left-hand side may be transformed to a double integral using (j20p . If we then 
exchange the order of integration and set C = cos^ a and S = sin^ (p in the evaluation 

/ ^ = d9 = K (VC + S-CS) 

Jo ^/{l-Ccos'^e){l- Ssin^e) ^ ' 

given in [58], we obtain an identity first derived by Glasser [38], namely 

K ^^1 - sin^ a cos^ c/)^ d(p = K(sin(a/2)) K(cos(a/2)) , (51) 
which was re-derived by Joyce and Zucker and recorded in [231 Eq. 2.3.5]. 

3.8 Sum rule 

Using the analysis of 03^0 above and in |10| we obtain a sum rule 
In Section 5 we conjecture an integral counterpart to this sum rule. 
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4 Four Bessel functions 

We may construct ^4 by folding one instance of 5*3, in (j22]l . with the discontinuity of 
another, to obtain 

POO 

Si{a,b,c,w) = / 2uD'^{a,b,u)Ss{u,c,w) du (53) 

Ja+b 

from which we obtain an evaluation of the on-shell value 



4arctanh(^^H_J , lylogjy) ^ vr^ 

Hi:=54(l,l,l,i)= / ^^^^^ ^" = 70 = ^ ^ 



2 



by the substitution u = y + l/y. 

4.1 The odd moments S4 2/C+1 

By differentiation of 



4 arctanh 



S4(l, 1, 1, ix) = / ^,{u+ir-x^y 



we evaluate 

1 , d d5'4(l,l,l,ia:) _ vr^ 



S4,3 = - + -;- 



X dx J dx 

and are then able to solve the recursion relation [26^ [60] for S4,2fc+i by the closed form 



S4,2fc+l = ^(^) h (55) 



with integers 

k 

J J V k-j J \ j 



j=0 

Integer sequence ([561) begins 

1, 4, 28, 256, 2716, 31504, 387136, 4951552, 65218204, 878536624, 12046924528 (57) 
and is recordecH as entry A2895 of the on-line version of [66] , which gives the recursion 

k\ - 2{2k - 1)(5A;2 - 5k + 2)bk-i + 64(A: - 1)^6^-2 = (58) 
and the generating function 

/o^2t) = j;6J- . (59) 
fc=o ^ ■ ^ 

We have verified that recursion ()58p reproduces the recursion for (j55p , which has the same 
form as for the odd moments in ([9]) . We note that in |36[ |43] the integers hk were related 
to enumeration of paths in three-dimensional diamond lattices. They also appear in a 



See http://www.research.att.coin/~njas/sequences/A002895 . 
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study \21l Eq. 6.10] of lattice magnetic walks. From the generating function Jq we see 
that they result from the constrained sum 

^-^ \plqlr\sl 

p+q+r+s=k 

The economical recursion in (1581) will be used in Section 5. 



4.2 Dispersion relation 

We adapt the dispersive method of [201 [2S1 ES] to two spacetime dimensions and take the 
discontinuity of (I53p across the cut with branch point w'^ = — (a + 6 + c)^, obtaining the 
dispersion relation 

S4{a,b,c,w) = / 5- 5 dv (60) 

Ja+b+c + 

with a discontinuity given by a complete elliptic integral K of the first kind, and hence 
by an arithmetic-geometric mean p] [25] , namely 



D4{a,b,c,d) = / 2uD2,{a,b,u)D^{u,c,d)du (61) 

J a+b 

(62) 



a+b 

C 

Q{a,b,c,d) I TT 



1' Q{a,b,c,~d) 



Q{a,b,c,d) AGM (Q {a, b, c, d), iVabcd 



where 

Q{a,b,c,d) := {a + b + c + d){a + b — c — d){a — b + c — d){a — b — c + d) 

is completely symmetric in its 4 arguments. 

In physical terms, D4{a,b,c,d) gives the volume of phase space for the decay of a 
particle of mass d > a + b + c into 3 particles with masses a, b and c in two spacetime 
dimensions. In 4 spacetime dimensions, one would obtain an incomplete elliptic integral 
for the area of the Dalitz plot, in the generic mass case; in two spacetime dimensions we 
obtain a simple arithmetic-geometric mean for the length of a Dalitz line [33j . 

Using this elliptic representation for 5*4(1,1,1,0) = V3(l,l,l) = 03^1 we obtain an 
evaluation of the integral 

^'z?(y)dy = C3,i = ^L_3(2) (63) 
with D{y) := 21)4(1, 1, 1, l/y)/y given by 



^ '^^IVITOK^J = ^ HeunG (-8, -2; 1, 1, 1, 1; 1 - 9,^) (64) 
V(l + 3y)(l-2/)3 2 ^ , , , , , , y ) K ) 

where the general Heun function, HeunG, satisfies Heun's differential equation as specified 
in Maple, see [6l]. Similarly, 

. o , 1 1 . 1 



I' D{y)y^dy = ^c^,3 = \L^3{2) 



6 
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In Section 5 we shall show that the HeunG representation of the elliptic integral ([61 
from quantum field theory, may be related to the hexagonal lattice sequence (j24p of crystal 
theory. For the present, we note that ([Mj) and (jUUj) provide the evaluation 



a D{y) ^ 7r2 



dy = S4,i = — . (65) 
1 - r 16 

4.3 The odd moments C4 2fc+i 

For V4 we have two representations. First, from the elementary arctanh function in (j22p 
we may compute 

V4(a, 6, c, (i) = / wS-i{a,h^w)S-i{c,d,w) dw (66) 

and easily evaluate 







4arctanh2 (^7^^) .i 4ylog2(y) ^ 7 



7o ^"(^ +4) Jo 1 - y 8 

by the substitution w = 1/y — y. Similarly, by differentiation of ([66]) . we obtain 

C4,3 = |C(3)-^. 

In general, all the odd moments C4^2fc+i are given by rational linear combinations of C(3) 
and unity, as shown in [TOl EH EQ]. Because of the mixing of C(3) with unity, we were 
unable to write a closed form for their rational coefficients in C4^2fc+i- 
The alternative folding, using ([62]) . is 

/>oo 

V4{a,b,c,d) = I wS4{a,b,c,w)S2{d,w) dw 

00 

2v Di{a, b, c, v)V2 {v, d) dw , 

I a+b+c 

which yields a novel formula for C(3), namely 

D{y)\og{y) 7 

2 ^ — = "^^.i = -C(3) (67) 











by the substitution v = 1/y. 



4.4 The even moment C4 

The analogue to (|66|) . for even moments, is 



00 



V4{a,b,c,d) = 2tt / S3{a,b,w)S3{c,d,w) dw (68) 







with a product of elliptic integrals coming from ()30p . Setting c = a and d = b, we obtain 
the intriguing case 



j^2 / /i«2+{a-6)2- 

°° i^o^(ai)i^o^(5i) dt = 2. 2^T?M^ (^9) 

13 



where the square of K may be replaced by a 3F2 series [15] for only part of the range of the 
integration, because of the restricted validity of the Clausen product (j36|) . In particular 
it was not at first clear how to evaluate the integral in 

C4,o = vr / K^{s{u6)de (70) 

obtained by setting a = 6 = 1 and w = 2tan0. 

The key to unlock this puzzle was provided by the trigonometric series 

00 

K(sine) = ^7„sin((4n + 1)6') (71) 

n=0 

with coefficients 

given by F.G. Tricomi [68] and recorded in [371 Section 13.8, Eqn. (8)]. The identity ([7T]) is 
valid for 7r/2 > ^ > 0. Thus, the integral in (I70p is easily transformed to a hypergeometric 
sum in 



2 00 4 °° ('2n^4 4 / 1 1 1 1 

•^4,0 4 7n 4 58n A 111 

n=0 n=0 



1 (73) 



4 4 4 V 1,1,1 

n=0 n=0 ^ 

by using the orthogonality relation 

/ sin((4m + 1)0) sin((4n + 1)0) dO = -5m,n ■ 
Jo 4 

It is instructive to split 04^0 = ^1 + -4.2 into the contributions 

Ai = 7^ K'^{smO)dO, A2 = TT K'^{cosO)dO, 
Jo Jo 

and use the Clausen product (j46p for the former. The result is ^1 = which proves 
that ^2/-4i = 3, as we had noticed numerically. A more direct derivation of this proven 
factor of 3 would probably be enlightening. 

We also note that Ai = |c4.o may be obtained by setting (f> = in the double series 



rw/A 3 °° /2n\3 / 1 1 _ 

vr/ K\sme)dO=^coscP^^2Fi( " cos' 

J4>/2 » „=o ^ V 2 

which is valid 7r/2 > (/> > 



4.5 Sum rule 

Combining (j73p with a more complicated sum in [10, Eq. 3-5], we get the discrete sum 
rule 

with one more central binomial coefficient and one more power of 7r/^/3 than dSSj). 
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4.6 The even moment C4 2 

Integer relation algorithms (see [SJ [71 [T^] for extended discussion) led us to conjecture 
that 

vr^ °° / 12 3 \ 

n=0 ^ V I ; / 

which we shall now prove. 

As before, we use Bessel's differential equation and operate on our master formula, in 
this case ([GS]) . before setting parameters to unity, obtaining 

r/Vcos4e-l 4cot^ + sin4^ d \ . m 



Then, using Tricomi's expansion (j7ip . we easily reduce 

rn/2 2 °o / i 

vr / K2(sin0)cos46'd6' = — 2 

Jo 16^0^ ^ + ^ 

to single sums of the form in (|75p . The term involving sin 46 gives a multiple of the same 
sum, using integration by parts. 

However, the term involving cot is more demanding. Integrating it by parts, we 
conclude that (I75p is equivalent to the evaluation 

7o sm20 4 ^ '"(n + l)2 

which we now prove by using a subtracted form for the differential of (j7ip . 
Defining (5„ := 4 — (4n + l)7n, we prove by induction that 

n—l oo 

"^6m = An-4n^-fn = l-^dm- (77) 

m=0 m=n 

Then from this we derive a subtracted series for 

dK(sin0) sin^6' ^ , ^ 

= >^5„ COS0-COS 4m + l 78 

at) cos ^-^^ 

by subtracting the trigonometric series [59] 

, x—^ sin((4m + 1)6*) vr , , ^ ^, , 

4 E L + i = 2 + ^ + ^^'^ ^^^^ 

m=0 

from (j7ip and then differentiating, to obtain 

dK(sin^) ^ ^ ^ 

— ^ _ sec e = - ) cos ( 4m + 1^ . 

m=l 

Adding cos Q to each side and using X^^^q 5m = l, as a consequence of ([77p . we obtain ([75]) . 
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By combining (j7T|) and (j78|) . we obtain the double series 

B 

^ - 70 = ^ (Mo,„ - Mm,„) 7„ 

m,n>0 

using the triangular array M with entries 

M„ n := - / cot e cos((4m + 1)61) sin((4n + 1)6*) dO 
^ Jo 

for non-negative integers m and n. These entries are very simple: Mm,n = 2, for m < n, 
Mm,n = 0, for m > n, and M„^„ = 1, resulting in 

^ oo / m \ oo 

^ = 70 + X] ( - 70 + 2 ^ 7n I = ^ 7nen 

m=0 V n=0 / n=0 

with the 7o terms cancelling, since X]m=o '^»" ~ ^- Here, 

oo 

en := + 2 ^ (5™ = (Sn^ + 4n + 1)7„ - 2(4n + 1) 

r?i=n 

appears by a change of the order of summation and is easily evaluated, by using ([77j) . 

Looking back to what needs to be proved, in ([76]) . we see that we now need to establish 
the vanishing of 

This is achieved by taking the N ^ oo limit of the explicit evaluation 

E7»(g^7„-.„)=2A.S„.„ = 0(1) (80) 

of a truncated sum, which is easily proven by induction. We note that the Gosper algo- 
rithm [62] in Maple failed to evaluate (jSOp . as written, since 6n and 8^ mix binomial and 
polynomial terms. If one separates these by hand then our compact result is verified. 

4.7 Further evaluation of integrals 

We also succeeded in separating B = Bi + B2 into the contributions 

Jo sm^ e 4 ^ n + 1 

n=0 

/n cos2 6' 4 ^ n + 1 

with a familiar factor of 3 multiplying the 4F3 series in the latter and a new constant from 
the singular value ki = sin(7r/4), namely 



n := ^Mi-3F2( 



r^(i) + i6r^(|) 

8 

The vr^ terms in Tl and the integrand of (I82p match, since J^^^ sec'^ 9 d9 = 1. 
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4.8 Further evaluations of sums 

Evaluations ([ST|) and ([52]) were discovered using PSLQ. To prove we may again use 
a Clausen product, with its first term subtracted. We then encounter the integrals 

2/ ^ — de = B{n--,-]+ 

70 sin^e V 2'2y 2n-l 

with n > 0. Here B is Euler's Beta function and yields yet another 4F3 series, which we 
eliminate by using the — > 00 limit of the summation 



N-l 



n=0 



1 



2n n + 1 



2iv3T^- = « + o - 



which was also proven by induction. To prove (j82p . we then subtract (I8ip from (I76p and 
use the N ^ 00 limit of the summation 



n=0 



— ^ + ^-^TT^ I = 16iV^7^ = 16 + O f ^ 
n + 1 (n + l)2y \^Af 



which was proven in the same manner. 

Thus one may undo the explicit evaluations of parts of (f82ll . in terms of tt and T values, 
and instead write 



^ ;-/^ K2(cos e) ^^_n'y^ {':f 2n + 1, vr^ An' + lOn + 5 

8 26" n + 1 32 ^ 28" in + 1)2 

n=0 72=0 







COS^ 



in terms of undigested hypergeometric series. 

More productively, we may use the explicit summations so far achieved to reduce 



16C4,o - 64C4,2 



— 4F3 



1111 
2 ' 2 ' 2 ' 2 
2,1,1 



(83) 



to a single 4F3 series. Comparing this last result with (|73p . we conclude that all moments 
C4,2fc can be expressed in terms vr and a pair of contiguous 4F3 series. An equivalent 
hypergeometric expression is 



C4,2 = ^ ^ 4 4F3 



1111 

2' 2' 2' 2 
1,1,1 



34F3 



1111 
2' 2' 2' 2 

2,1,1 



37r2 
16 



(84) 



For comparison we repeat 



TT 



C4,0 



4-f 3 



1111 
2 ' 2 ' 2 ' 2 
1,1,1 



(85) 



4.9 Relation to Meijer's G-function 

A generalization of (j69p is derivable in terms of the Meijer- G function \63\ Vol. 3] . For 
example, we have 



I{a,b,k):=l t''K^{at)K^{bt) dt 
/o 



vr 



1-fc 1-fc 1-fc 1 
2 2 2 2 



(86) 
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which may be proven by making two copies of the representation 



2 
H 
2 



2 



and integrating them with weight t to obtain 
I{a,b,fi + u + 1) 



IT 



13 



ah'' 



2 

H 

2 



13 



2 



dt. 



Then we use Meijer's result for the integral of the product of two G-functions to obtain 

62 



I{a,b,^ + u + l) 



^33 

8^^+25.-^44 



2 



2 



E 
2 



u+1 
2 

2 



The apparent freedom in the choice of parameters /x and i' is demystified by formula 
8.2.1.15 in [631 Vol. 3], which shows that multiplication by a power of the argument of a 
G-function is equivalent to adding a constant to all its parameters, as in the example ()38p . 
We resolve this redundancy by setting fj, = k — 1 and 1/ = and hence prove (j86p . Similarly, 
we obtained 







t''Ko{at)Kl{ht)dt ^' 



a 



k+l 



^32 ^462 



1-fc 1-fc 

2 2 





1 

2 





(87) 



In [To] the special cases of (i86]l and (fHTI) with a = 6 = 1 were studied numerically, using 
Adamchik's algorithm [2]. This algorithm converges quickly in the case of c^^k-, obtained 
from (I87p with an argument of 4. But in Maple it is painfully slow in the case of 04^^, 
obtained from (j86p with unit argument. Numerical evaluation of our new hyper geometric 
results (f84|l and (f85|) goes far faster, with Maple. 



4.10 Another continued fraction 

Again we may derive from pi]] that 



C4,2 
C4,0 



efll 



(27V - 1)^ 



where e(A^) := A^(20A^2 _)_ 3) ^^d the ratio 8c4^2/c4,o niay be made explicit from (j84|85p . 



4.11 The even moment S4 







The odd moment 84^1 relates directly to quantum field theory; it is the two-loop on-shell 
equal- mass sunshine diagram in two spacetime dimensions. No such meaning attaches to 
S4^oi it is hard to think of a physical application for this moment. However, we found a 
rather pretty formula for it, which we record as 



S4,o := / h{t)K^{t)dt 
10 



n/2 



K(sin6l)K(cos6l)d6'. 
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This amusing twist of the integral (j70p for 04^0 follows from Nicholson's integral represen- 
tation (721 13.72, Eq. 3] of the product 

2 r/2 



2 r' 

Io{t)Ko(t) = - Ko{2tsma)da. 
71" Jo 



Substituting c = 2 sin a in (|Mj) and using the appropriate reduction (j36l) to a product of 
K values, we obtain 



/ Io(t)K^(t)dt= / ^ ^ da= / K(sin(a/2))K(cos(a/2))da. 
Jo Jo 2sma Jq 

Setting a = 26, we prove ([89|) . Then Pari-GP gives, in a tenth of a second, 64 digits of 

S4,o = 6.997563016680632359556757826853096005697754284353362908336255807 . . . 
A corresponding twist of the sum (j73p for C4 comes from Tricomi's expansion (j7ip and 



7r/2 

sin((4m + 1)61) cos((4n + 1)6^) d9 



4m + 4n + 2 
This yields the double sum 



E 7m7n _ 2 1 ^ 

4m + 4n + 2 " ^« ^ 



m,n>0 n=0 

with 7n in ([72]) and 



A„ := A y = 1 + y 



7„ 4m + 4n + 2 16m^ — 1 

m=0 m=l 

where the latter form was proven by Maple. Hence we obtain 

n=0 \ fe=0 / 

By way of comparison, we note the simpler evaluation 

/ij' — — 

from the closed form for C2,2fc in © and the expansion 

4.12 Tabular summary 

In Table 1 we recapitulate the key discoveries for the moments c^^k •= fQ^t^^oi't)^^ 
with n = 3, 4. The results for the even moments 03^2^ and 04^2*: are new. The table may 
extended by using the recursions dH]) and ([9]). 
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C3,0 
C3,l 

C3,2 
C3,3 



3r-(i) 

327r22/3 
967r22/3 916(1) 



3J^ 2 



/111 


9 


1 2' 2' 2 







288 



3-r 2 



111 
2' 2' 2 

2,2 



C4,0 
C4,l 
C4,2 
C4,3 



2n\4 



4 oo I2n\ 
n=0 

k(3) 



TT 



4-^3 



1111 
2' 2' 2' 2 
1,1,1 



TT 



1111 
2' 2' 2' 2 



16 ''''V l.l-l 

32^^ ^ 16 



64 



4-r 3 



1111 
2' 2' 2' 2 

2,1,1 



37r2 
16 



Table 1: Evaluations for Cn,k with n = 3,4. 

5 Five Bessel functions 

Little is known for certain about integrals involving 5 Bessel functions. However, there 
are some remarkable conjectures arising from studies in quantum field theory [55^ I56j. 



5.1 Conjectural evaluations of Feynman diagrams 

In [55], Stefano Laporta developed an impressive technique for numerical evaluation of the 
coefficients of the Laurent expansion in e of Feynman diagrams in D = 4 — 2e spacetime 
dimensions. Here we are concerned with just one of the many diagrams that he considered, 
namely the dimensionally regularized 3-loop sunrise diagram with 4 internal lines: 




d^Pi d>2 d^p3 



N{pi)N{p2)N{p3)N{q -PI-P2- Ps) 



q-q=w^ 



where N[p) := p ■ p + 1 is the inverse propagator of a scalar particle with unit mass and 
momentum p. In the on-shell case, the Laurent expansion found by Laporta has the form 

S,{-lA-2e) _ 2 22 577 

with a numerical value ^205 ~ 21.92956264368, for the finite part, given in equation (205) 
of [55]. Subsequently, in equation (21) of [56], this constant was conjecturally related to 
products of elliptic integrals of the first and second kind, with a numerical check to 1200 
decimal places. 
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In a entitled "Reciprocal PSLQ and the Tiny Nome of Bologna", Broadhurst 

observed that Laporta's conjecture may be written rather intriguingly using the constant 

C := -- ( 1 - ^ ) ( 1 + 2 > ' exp ( -n^vrVls) ) (93) 



and its reciprocal 1/C, in terms of which he wrote Laporta's conjecture as 
Broadhurst further conjectured that the odd moments 

l-OO 

Jo 

are linear combinations of vr^C and ir'^/C, with rational coefficients, and in particular that 

(95) 

13C-J^) (%) 

43C-i^) (97) 

with higher moments obtained by a recursion of the form (jlOp with polynomials p^^i given 
in (fTTIl . We have checked these 3 conjectures to 1200 decimal places. 
In the course of this work, we discovered that the moments 



•55,1 








•55,3 








•55,5 


m ( 


7r2 





Jo 



follow an uncannily similar pattern. If and are the rational numbers that give 
•S5,2A:+i/'^^7 conjecturally, as q^C — rf^/C, then we found that q^C + r^/C gives the value 
of 2t5^2fc+i/(V^vr). We checked 60,000 decimal places of the resultant evaluations 

(98) 

1 



2*5,1 




C 








2*5,3 




(- 








2*5,5 




(- 


Vl57r 




1 15 



10C7 
40C7 J 



(99) 



(100) 



for which we eventually found a proof, presented in subsection 15.101 

Finally, by doubling one of the masses in the Feynman diagram corresponding to ts^i, 
we arrived at the conjectural evaluation 

?r«i 1 

tli{t)Kl{t)Koi2t)dt =^ J^K^K'^ (101) 
which has been checked to 1200 decimal places. 



•^Zentrum fur interdisziplinare Forschung in Bielefeld, 14th of June, 2007. Displays available from 
[http : / / www ■ physik . uni-bielef eld . de/igs/schools/ZiF2007/Broadhurst . pdf ^ leading to 200,000 decimal 



places for C5,i and 05,3 in http://paftp.open.ac.Uk/pub/staffjftp/dbroadh.u/newconst/V5AB.txt 
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Notation. In the 8 evaluations ([M|) to (jlOip we have used the device = or = 
to distinguish the cases that remain unproven from the 3 cases in (j98p to (jlOOp . which 
we were eventuahy able to prove. Some of the labels n = 1 ... 8 will recur, as we give 
equivalent forms of these conjectured or proven evaluations. 



5.2 The odd moments ^5,2^+1 

Evaluations ()98|99|100p were easy to check to high precision, thanks to our closed form ([55 
for the odd moments 54^2^+1- By expanding one of the functions 



n=0 ^ ' 



in the integrand, t^'^^^/g (0^0(^)5 °f ^5,2fe+i5 we obtain a rapidly converging sum in 



n=k 



W = 4-V5:6„(^,,^) (102) 



in terms of the diamond lattice integers ([56]) . To relate t^^i to a product of complete 
elliptic integrals we use Jacobi's identity 



^ = ^3(.) := f; (103) 



with a nome related to k hy q = exp(— 7rK'(fc)/K(A;)). Specializing to the singular 
value [21] 

, (2-^/3)(^/5-^/3)(3-^/5) 
= ^2 

with the "tiny nome" gi5 := exp{-^T^/T5) « 0.000005197, we obtain from ([93D 

with i^is := K(A;i5) and K^j^ := ^.{k^/^), where 

= ^ 

yields the larger nome ^5/3 := exp(— Try 5/3) = q-^^ . Thus evaluation ([98|) amounts to 

«=^.=^f:^ = ^^wvv= (104) 



n=0 

with a summand 6„/64" = 0(n^^/^/4"), from (I56p . giving rapid convergence. By taking 
10^ terms, we checked ([98]) to 60,000 decimal places, using the recursion ([58]) for the 
diamond lattice sequence bn- Our closed form in ([104p resulted from paying diligent 
attention to a footnote in [561 P- 121], which led us, eventually, to discover and prove the 
connection between quantum field theory and these diamond lattice integers. 
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We remark that our evaluations of S5^2fc+i and ^5,2^+1 naay be expressed in terms of 
r values, using the corresponding evaluation of Kf^ in [7], as was remarked by Laporta 
after Broadhurst's talk (see footnote 3). For example, we may re- write the conjectural 
evaluation for S5 3 in ([96]) as 



^ r^^mimt a "r(A)r(A)r(A)ra) _ ra)r(li)r(S)r(M) 

2 Jq 30 15 

(105) 

which contains all 8 values of r(n/15) with n £ [1,14] and coprime to 15. Then the 
counterpart for t^^s in (|99|) may be written as 

- ^[61 i3r(^)r(A)r(A)r(A) r(X)r(ii)r(i|)r(i|) 

4^/3 ^ 64" 303 ^ 15 y^""^) 

^ n=l 

by the remarkable sign change discovered in our present work and the relation to diamond 
lattice numbers in (I102p . 

5.3 Double integrals 

The moments cs^i, ss^i and t^^i are easily expressible as double integrals of elementary 
functions. 

For the 4-loop vacuum j56j diagram 



00 



V5{a,b,c,d,e) := / t Ko{at)Ko{bt)Ko{ct)Ko{dt)Ko{et) dt 
Jo 

in two spacetime dimensions, we obtain the double integral 



POO /'OO 

V5{a,b,c,d,e) = / xyS3{a,b,x)D3{x,y,ic)S3{d,e,y) dx dy 
Jo Jo 

by grouping the internal lines with masses a and b to give a total momentum with norm 
and those with masses d and e to give a total momentum with norm y^. Then the 
coupling term 

1 r de _ 1 - D {x y ic) 

IT Jq x2 + 2xy cos 6* + + y^(x + y)^ + c2y^(x - y)2 + 

comes from an angular average in two Euclidean dimensions. Setting the five masses to 
unity, we obtain 

_ /■°° r°° 4arcsinh(x/2) arcsinh(y/2) dxdy 

""'''"70 Jo V(4 + 2;2)(4 + y2)(l + (2; + y)2)(l + (a; _ y)2) 



where we have converted the arctanh function of (j22p to a more convenient arcsinh 
function, in the equal-mass case. 

Similarly, the 3-loop sunrise diagram 



00 



S5{a,b,c,d,z) := / tKo{at)Ko{bt)Ko{ct)Ko{dt)Jo{zt)dt 
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yields the double integral 

S5{a,b,c,d,z) = / xyS3{a,b,x)D3{x,y,ic)D3{y,z,id)dxdy 
Jo Jo 

by cutting one line in V5 and setting the norm of its Euclidean momentum to z'^. Setting 
the 4 masses to unity and analytically continuing to the on-shell point = — 1, we obtain 

/■°° 1"^ 2arcsinh(x/2)dxdy ?[2[ ^ K 

~ Jo Jo V(4 + x2)(4 + + + y)2)(i + _ y)2) — ^Tff 15 5/3 

(108) 

whose conjectural evaluation is given by ([95|) . 
We then define 

POO 

T^{w,a,h,c,z) := I tJ(j{wt)Ko{at)KQ{bt)Ko{ct)Jo{zt)dt 
Jo 

as the angular average of the diagram obtained by cutting two lines in V5 and setting the 
norms of their momenta to w'^ and z^. Hence we obtain the double integral 

Tr,{w,a,b,c,z) = / xyD3{w,x,ia)D3{x,y,ib)D3{y,z,ic)dxdy 
Jo Jo 

which leads to 

, _ r r dxdy ^ K nnq^ 

/ / , ou. , , , , " 7^15^5/3 (109) 



JO 



V(4 + x2)(4 + y2)(l + + y)2)(l + (x - y)2) 



whose evaluation is given by ([U5|) . If we multiply pU9|) by 4, we recover the double- integral 
discovery reported by Laporta in Eqs. 17, 18 and 19 of [56j. 
Finally, by doubling one of the internal masses, we obtain 



00 POO /"OO 

tliit)Klit)Koi2t)dt = 
Jo 



/o V(4 + x2)(4 + y2)(4 + (x + y)2)(4 + (x - y)^) 
and hence the conjectural evaluation 

r a i (110) 

A Jo + l-'-'Xl + + (l- + !/)-)(! + (l- - 3 ^ " 

after rescaling x and y by a factor of 2 and invoking (jlOip . 

Evaluation (jllOp resonates with the proven evaluations and (H5]) in Section 3, 
where we found that removing (1 + (x — y)^) from the square root in (jllOp doubles the 
value of the integral and that removing (1 + (x + t/)^) multiplies it by 4. 

We were unable find a transformation of variables for the double integrals in (jl09p 
and (IllOp that suggested their evaluations as products at the singular values ki^ and k^, 
respectively. In the next 3 subsections we show how to express ()109p as a single integral 
of a complete elliptic integral, in 3 rather different ways. 

As observed in a footnote in [561 P- 120], entry 3.1.5.16 in [631 Vol. 1] is intriguing: for 
real parameters, ki and /c2 with k\ + k"^ < ^ one has 

/-/2 _2K(a)K'(/3) 



./o 



1 - kl sin2 6* - A;2 sin^ <p 1 + \/l - ^: 
24 
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where 

a := ; and := , -. 

Perhaps there are imphcations for (|110p from this general form. As discussed in |19^ [39] 
one can estabhsh a more recondite counterpart for 

rTT/2 pTT/2 



5.4 Single integrals from polar coordinates 

We may recast evaluations (|lU9p and piU|) by transforming the more general double 
integral 

dx dy 



oo roo 



/o Jo V(c + ^')(c + y2)(l + {x + y)2)(l + (x - y)2) 

to polar coordinates. With x = r cos and y = r sin 0, the product of the first two factors 
in the square root gives a term linear in w := cos 40, as does the remaining product. An 
angular integral of the form 

dw ^ K ^ VA^ + B^ - l\ 



^{1 - w^){2A'^ - 1 - u;)(252 -1 + w) AB y AB J 

results, with A = 1 + 2c/r'^ and B = 1 + l/r^. Transforming to z = j (1 + r'^), we obtain 
the single integral 



J(.)=/'k i^^i + MZpIElZl) (112) 







z + 2c(l -z) ) z + 2c(l - z 

Setting c = 4, we transform evaluation (I109p to 



Y^(13- 12z)(5 -4z) \ dz ^ 1 







K ^ -Ki^K./^. (113) 

> 8 - 7z 8-7z 4 i& 5/3 1. ; 



Setting c = 1, we transform conjecture (jllOp to 







/ zV5-4z\ dz ?[8] 1,,., ,,,,, 
K \ ^ — — = -i^3i^i/3 (114) 



since, by definition of a singular value, i^i/s = -fC3. A hyper geometric version of (|114p 
may be obtained by writing its left-hand side as the ornate triple sum 



/ 2n + l,2n + l 
n^({2n)iy^^^, j''[ 2n + 2 + m _ 
4^V22"n!/ ^ (n - m)\(2n + 1 + m)l ^ \ 1,1 

71=0 m=0 \ / \ / V \ 



1 I 



2 ' 7r«l 7r2 /111 
3F2 ( 2' 2' 2 



We computed the single integrals in (jll3p and (jll4p using Pari-GP, which provides 
an efficient agm procedure, for evaluating the complete elliptic integral K, and an efficient 
intnum procedure, for the numerical quadrature. In each case, we confirmed the evaluation 
to 1200 decimal places but were none the wiser as to its origin. 
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5.5 Single integrals over discontinuities 

Seeking illumination, we turned to integrals over the elliptic integral (j62|) in the disconti- 
nuity coming from the Dalitz-plot integration in (j6ip . 

For the 4- loop vacuum diagram V5 we may fold D4 from (j62p with V3 from (|28l) . to 
obtain 

I" 00 

V5(o, 6, c, d, e) = / 2?;D4(a, 6, c, t')V3(t>, d, e) dt; . (115) 

Ja+b+c 

With unit masses and y = 1/v, this gives the moments 
D{y)L{y) ^ 

C5,i = / , ^ dy (116) 
Jo VI - 4y2 

with a complete elliptic integral in D(y), from (j64p . and a dilogarithmic function 



m Il, - (118) 



2 \ v^l - 4y2 + 1 / 2 \ ^l-4y2-l 



7r2 



with (|118p coming from (|28p and the reduction to a single convenient Li2 value in (|119p 
obtained by transformations in ^57^ A. 2.1] and noted in ^34^ Eq. 3.21]. Further differenti- 
ations yield an even lengthier integrand for 05^5. However, we expected the odd moments 
C5,2fc+i with /c > 1 to be expressible as Q-linear combinations of 05^3, cs^i and unity. Hence 
we lazily used PSLQ to arrive at our ninth conjecture 

?[9] 76 16 8 

C5,5 = Y^C5,3--C5,l + - (120) 
with a final rational term originating, presumably, from the analytically trivial input 

r fKomtit) dt = hm = i (121) 

Jo 5 5 

to the richer (and more challenging) recursions considered in the talk of footnote 3, which 
dealt with integrals of products of powers of t, Ko{t) and the derivative -fCo(^) ~ —Ki{t). 
From (I116p . (I117p and the eminently believable conjecture (I120p . higher odd moments are 
obtainable from (jlOp and (jlip . Unfortunately we lack a more explicit evaluation of the 
single integrals for 05^1 and 05^3. 

For the 3-loop sunrise diagram, the corresponding folding is 

POO 

S5{a,b,c,d,w) = / 2v D4{a,b,c,v)S3{v,d,w) dv (122) 

Ja+b+c 

from which we obtain the on-shell value 



S5 1 = [ — arctanh | -1 / ^| dv = — ^i^isi^t;/^ (123) 

' Jo yr^^ vV 1 + 22/7 2VT5 ^ ' 
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and, by use of Bessel's equation, the more complicated integral 



= i (^w [ VI -V - ' ^ ^ ' ''-'^ 

for the next odd moment, conjecturally given by (|96p . Then PSLQ gives 

?[2,3,4] 76 16 
S5,5 = 77 55,3 - 77 S5,l (125) 
15 45 

which presumably results from a partial integration in the even richer (and even more 
challenging) recursions for integrals of products of powers of t, Io{t), Ii{t), Koit) and 
Ki{t). We used Pari-GP to evaluate the dispersive single integrals for the first 3 odd mo- 
ments S5,2fc+i) with A; = 0, 1, 2, to 1200 decimal places, and hence checked the conjectured 
evaluations (j95p . (j96p and (j97p and their consequent integer relation (jl25p . to this high 
precision. 

Similarly, the folding 

T^{u,a,b,c,w) = / 2v D4{a,b,c,v)D^{u,w,w) dv (126) 

Ja+b+c 

gives, with y = 1/v, the on-shell, unit-mass result 

= / /, , = -Ki5K^/3 (127) 

Jo V 1 -"^y ^ 

and, by use of Bessel's equation, 

_ f UyHl-2y'+4y^)D{y) 
''''-Jo (1 - V)^/2 ^'''^ 

for the next odd moment, given by (|99p . Then PSLQ gives 

^5,5 — 77*5,3 — 77*5,1 (129) 
15 45 

which was checked to 60,000 decimal places, in the more convenient sum 

y ( lQn\n -If-^-^ 4n2 + ^1^0 (130) 
^ 64" V 15 45 / ^ ' 

n=0 ^ ^ 

over the diamond lattice integers in ()102p . 

By setting u = 'w = \^a = b = \ and c = 2 in (I126p and transforming to y = 4/f, we 
obtain the representation 



/ ^ i/y +v ^ ^ -i^3i^i/3 (131) 

70 V(2-yP(2 + y)(i + y) ^ 3 ' 

as the dispersive counterpart to (jll4p . but still have no better idea as to why the third 
singular value occurs. 
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5.6 Relations between elliptic integrals 

We found a third way of stating our conjectures in terms of integrals over K. This arose 
from the double-integral representation 



2 arcsinh(x/2) dx dy ?[2] tt 
Jo ^{4 + x^){i + y^)Rix,y) ~ 2^15 



«5,i=/ / ... . .... : . = :r^i^i5i^5/3 (132) 



with the product (1 + (x + y)^)(l + (x — y)^) in (llOSp replaced by the non-factorizable 
term 

■ (x + i)2 + 4 + y2 



R{x,y) 



3 + x2 + y2x2 



We obtained (jl32p by introducing a Feynman parameter, z, to combine two unit-mass 
propagators, with momenta q — "p and r — in the integral 

1 dz 



AB Jo {zA + {1- z)By 

with A := A^(g — p) and B := A^(r — q). Here, p is the combined momentum of the 
other two internal lines and r is the external momentum of the 3-loop sunrise diagram. 
Integration over the two-dimensional vector q leaves an integral over z with a denominator 
z{l — z){r — p) ■ [r — p) + 1 that is symmetric about the mid-point z = ^. Integrating over 
the angle between p and r, setting p ■ p = x'^ and analytically continuing to the on-shell 
point r ■ r = —1, we obtain ()132l) by making the transformation z{l — z) = 1/(4 -|- y^) 
which maps z = to y = oo and the mid-point z = ^ to y = 0. 

This method provided an analytical advance, since it proved the equality Ei{u) = 
E2{u) between the elliptic integrals 

Ei{u) := - / (133) 

2 Vr;(4 + v){l + 2u + 2v + {u - vf) 

E2{u) := / , (134) 
Jo y^v{4 + v)i4u + {3 + u + v)^) 

where we have transformed to u = x^ and v = y^. When Maple was asked to evaluate Ei 
and E2 it printed different expressions, each containing an incomplete elliptic integral 

1 

F{sm(P,k):= =de 
Jo V 1 - A;2 sin^ 9 

with a relation between k and (p which we reduced to the form 

(1 — cos (p) sm (p 

The incomplete integral was eliminated by computing {2Ei + E2)/3, for which Maple 
gave an expression involving only the complete elliptic integral K(/c). By this means, we 
proved that 

F(sin0,A;) = |k(A:) (136) 
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whenever k and are related by condition (jl35p . Moreover, by using [1., 17.4.13], we 
proved that this condition also gives the evaluation 

F{l-cos(j),k) = ^K{k) (137) 

which we shall use in the next subsection. 
Our result for E{x) := £'i(x^) = -^2(2;^) is 



with 

3arctan(x) — arctan(x/3) 
a(x) = . (139) 

This then gives 

2E{x) . , /X\ _ ?[2] TT 

as the non-dispersive counterparts to the integrals ([123]) and (jl27p over the different 
complete elliptic integral in D{y), with a closely related HeunG function in (j64p . 

We note that -E(x) contains a factor of 1/3 from evaluating {2E1+E2) /3. This ensures 
that we correctly reproduce 

S3,l = E{0) = ^ = ^-3(1) 

using a{x) = Ax/3 + O(x^). The relationship 

gives -E'(x) = log(x^)/x^ + 0(log(x)/x^), as expected at large momentum. 



5.7 Expansions near singularities 

The HeunG forms for D{y) and E{x), in (j64p and (|138p . yield expansions near the regular 
points y = \ and x = 0, respectively. However, these regular expansions were not needed 
in our numerical integrations, since the agm procedure of Pari-GP is highly efficient for 
the evaluation of K(A;) when k is not close to the singular point at A; = 1. What we really 
need, and eventually found, are the expansions 

00 

^(/ifc + afclog(y2))y2fc (143) 

k=0 

^ hk-ak log(x2) 

(-x2)fc ^ ' 



Djy) 

x^£^(x) = - 
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that isolate the logarithmic singularities as y ^ and x oo, respectively. Here au is 
the hexagonal lattice integer of (j24p and is determined by the differential equation for 
D{y), or equivalently E{x), which yields the recursion 

{k+lfhk+i-{Wk'^ + Wk+3)hk+9k'^hk-i = -2{k+l)ak+i + W{2k+l)ak-l8kak-i (145) 

with a starting value /iq = 0. We note that (jl43p converges for < y < ^ and (1144 1) 
converges for x > 3. We were alerted to role of the hexagonal lattice integers by the 
regular expansion 



vr v-^ / -x^ 



which is valid for |x| < 1, since = 0{9^/k) for large k. 

We were able to solve the recursion (|145p in closed form, by considering the moment 

T4{u,a,b,v) := / tjQ{ut)KQ{at)KQ{bt)jQ{vt)dt= / 2w D3{a,b,'w)D3{u,v,iw) dw 

Jo Ja+b 

which yields, in general, an incomplete elliptic integral 

r°° dr 
T4{u,a,b,v)= , (147) 

J(a+6)2 ^{x - (a + bY){x - (a - bY){x + (u + vf){x + (u - vf) 

by transforming to x = w^. In the special case with a = 1, 6 = \jy > 1 and n = u = i, 
Maple gave the evaluation 

. , 2y^F(l-cos(/>,fe) 
TA\, 1, l/y, i) = — , 

V(i + 3y)(i-y)^ 

with 



A. ( '^y \ V. (i-3j/)(i+j/)3 

(/) = arccos , k = \ — — — . (148) 

\l + yj V (l + 32/)(l-2/)3 ^ ^ 

This relation between k and </> is precisely as in ()135p . Hence, using (jl37p . we obtain 

=r4(i,l,l/y,i) := r tll{t)KQ{t)Ko{t/y)dt (149) 

since D{y) in (164p contains the complete elliptic integral K(A;) with k given in p48p . Then 
we used the expansion of /o(0> ([92|) . and of [1, 9.6.13] 

/ j-fc \ 2 A: 

i^o(t) = -(log(t/2) + 7)/o(t) + E 2] n ' 

fc=l ^ n=l 

where 7 is Euler's constant, to obtain from (jl43p and (jl49p the closed form 

^.^EG)^a)E^ 

j=l ^"'^ \ ^ n=0 

which is a harmonic twist of the closed form (j24p for the hexagonal lattice integers a^. 
Hence we obtain an integer sequence for '■= klh^/i, with A; > 0, beginning with 

1, 13, 263, 7518, 280074, 12895572, 707902740, 45152821872, 3282497058384... 

Thanks to this sequence, we were able to evaluate all the single integrals over D{y) and 
E(x) in this work to 1200 decimal places, since we had very good control of logarithmic 
singularities near the endpoints, as y ^ and x — > 00, respectively. 
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5.8 A modular identity from quantum field theory 

A careful analysis ^33j of the Dalitz plot shows that the relation (|135|) between k and (j) 
implies that [33l Eq. 5.17] 

Z(F(sin (f), k),k) = \k'^ sin^ (p (152) 

where Z is Maple's JacobiZeta function. Combining our new finding ()136p with this, we 
obtain from [H 17.4.38] an elementary evaluation of 

^ V ( ' ^ ] = -t-^ - -t-' (153) 

4K(A:) ^V(9~^""^-9^'"+^) (g-3«-2 _ g3n+2) y 2 6 ^ ^ 

where q is the nome associated with k and t = tan((/>/2) is determined algebraically by 

i / y 2(l-2k^)r 



(4/0^(1 - /c^)) 3 , s = 2A:\/r + A:2 + 2A;^ + 



+ r + r 



2 



^ _ / k-Vr + k^ + ^ y 
\3k + VfTW-VsJ 

In particular, for k = 1/V2 we obtain t = {2/\/3 — 1)^/^. 

A modular setting for this result is provided by Jacobi's identity 

OO 2 

Vkes{q)=e2{q):= ^^"""^^ (^^4) 

n.=— OO 

with ^3(9) related to K(fc) by (jlOSp . Summation of the Lambert series in (jl53p gives 
02{l^^'^) /G2{q^^'^), which is a result known to Ramanujan. Thus we have proven the 
modular identity 

where c := t ^ is the unique positive root of the polynomial 

S{x) = 3 + 8 (^1 - X + 6x2 - ^4 ^^5g) 

and indeed 1 < c < 3. The real roots are —0'^iq^^^)/O'^{q) and 30|(g'^)/0|(g), as has been 
known since Joubert and Cayley [Ml (4.6.14)]. Hence (I155p devolves to 



Such modular identities are machine-provable — in principle and in practice — by com- 
puting that sufficiently many terms of the Taylor series agree. This is the so-called "modu- 
lar machine" |25[ §3]. In this case confirming 1000-term agreement is more than adequate 
to prove (jl57p : as takes seconds in Maple. A conventional proof can be pieced together 
from [24, Thm. 4.11]. Our proof came from quite another source: the Lorentz invariance 
of quantum field theory in two spacetime dimensions, which enabled us to prove that ()133p 
and (I134p are identical. 
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5.9 A discrete sum rule 

It seemed reasonable to try to prove the evaluation for t^^i in one of its forms 

_ H Diy)dy m I ^ ^ 



/ E{x) dx A5] 1 ^ .^.„^ 
^M=( = 4i^i5i^5/3 (159) 



with D{y) given by dMD and E{x) by (fT38]l . 

For example, by using [23l (2.35), p. 119] one may deduce that 

/ / = - Ki^K^,^ 160) 

io Jo ^64 - (16 - sin^ 0) sin^ 8 ^ ^ 

and then hope to relate some such recognizable double integral to the single integrals 
in (fT58|) and (fT59|) . 

In fact, the integral forms (|158p and (jl59p resisted prolonged and intense efforts to 
find a proof. A break-through came from the sum in (jl04p . which we rewrite as 

= 16^ ^^^^^^ 

k=l 

whose right-hand side we shall relate to the first of the generalized Watson integrals |49l [50] 

.TT r-K d6lid6l2d6l3 , , 



with 



n n\ COS 6*1 COS 6*2 + COS 6*2 COS 6*3 + COS 6*3 COS 6*1 

wi(6'i,6'2,6i3) = (163) 



in the case of a face-centred cubic (f.c.c.) lattice 



(Q a Q\ cos 6*1 cos 6*2 + cos 6*3 
W2(6'l,fc'2,fc'3) = ^ (164) 



in the case of a simple cubic (s.c.) lattice 

1^3(91, 02,03) = cos ^1 cos 6*2 cos ^3 (165) 
in the case of a body-centred cubic (b.c.c.) lattice [17| and 

„ „ . 1 cos 6*1 cos 6*2 + cos 6*2 cos 6*3 + cos ^3 cos ^1 

W^4(fc'l,6'2,^3) = ^ (166) 



in the case of a diamond lattice 

In 1971, Joyce gave the notable f.c.c. evaluation [56', Eq. 4] 



12 K(fc+(z))K(fc_(z)) 
H^i(^) = ^ ^Tl ^16^) 



with 



2,^^_l, 2^/3 z ^3 (3 - z)(l - z)^/2 
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Our progress resulted from the intriguing observation that at z = ^ this gives fc_(l/5) = 
ki5 and A:+(l/5) = k^/^. Hence we obtain 

which reduces (|16ip to the discrete sum rule 



64fc 60^ 

k=l k=l 



where the f.c.c. lattice integers 

pTT 

fk'=—^ / / (cos ^1 COS 6*2 + COS 6*2 COS 03 + COS 6*3 COS 01 ) d6'id6'2d6'3 (171) 

TT-^ Jo Jo Jo 

give the Taylor coefficients of the expansion 

00 ^ 

Wi{z) = Y.fk{{^) (172) 

fc=0 

for \z\ < 1. In the next subsection, we show that (|17U|) derives from a more general set of 
sum rules. Then, in subsection 15.111 we shall expose a cubic modular equation, implicit 
in (fTHHD . 

5.10 Sum rules for diamond and cubic lattice integers 

The first few terms in the sequence for the f.c.c. lattice integers (I17ip are 

1, 0, 12, 48, 540, 4320, 42240, 403200, 4038300, 40958400, 423550512, 4434978240 

(173) 

for n = 0. . . 11. The values up to /s = 4038300 are recorded in [70l Table 1] and the 
next integer, fg = 40958400, was given by Domb [36] and recordeco in entry A2899 of 
the on-line version of [66j . which provided us with no closed formula for fk- By way of 
contrast, the s.c, b.c.c. and diamond lattice expansions 



k=o ^ ^ 

= E J (s) ("5) 



fc=0 
00 

W,{z) = Y.bk{{^y (176) 

k=0 

lead to explicit expressions for the integer sequences in entries A2896, A2897 and A2895, 
respectively, of the on-line version of [66]. As in the f.c.c. case p72p . expansions ()174p 
to (jl76p are valid for \z\ < 1. For convenience, we repeat here the closed forms 



See http://www.research.att.coin/~njas/sequences/A002899 . 
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for the hexagonal and diamond lattice integers, previously given in (j24p and (j56p . The 
hexagonal lattice integers appear in the s.c. lattice integers {^^)ak of expansion (|174p . 
We also note the terminating hyper geometric series 



for the diamond lattice integers. Likewise 

1 

2' 

1,1 

We were able to derive the closed form 

k 



= E ( J (178) 



j=0 



for the f.c.c. lattice integers, by noting the similarity between (1163P and ()166p . which leads 
to the functional relationship 

^ ^4 f TT^l =TrrWl (t?) (179) 



4iV + 4 \N + 1J m \N 

between the Green functions for the diamond and f.c.c. lattices. Then for > 3 we may 
expand each side, to obtain the sum rule 

^ (iivTI)^ " ^ (4iV)^ 

and derive the closed form ()178p for the f.c.c. lattice integers defined in (I17ip by further 
expanding the left-hand side in powers of 

Setting A'" = 3 in (I180p . we prove the sum rule 

oo , oo J. ^ ^ 

ETii = E^ = l^»'.W = i^^^V.V3 (181) 

k=l k=l 

from Watson's evaluation [71] of Wi(l), at the singular value k^. 
Setting = 15 in (jlSOp . we prove the sum rule 

i:|^-E^-i^'^'i(l)-i^'^.^%3 (182) 

k=l k=l ^ ^ 

at the singular value ^15. Hence we have proven (|170p and all the other forms of our initial 
conjecture (I98p for the moment ^5^1. By taking up to 4 differentials of Wi(3/A^), before 
setting A'' = 15, one may also prove the evaluations of ts^s in ([99]) and ts^s in (jlOOp . using 
the evaluations of the elliptic integrals E(A;i5) and E(A;5/3) given explicitly in the talk of 
footnote 3 and derivable from identities in [24]. 

These eventual proofs of our conjectures for the moments ^5,2^+1 came from noting a 
parenthetical remark in [48 » p. 601] to the effect that the Green function for the diamond 
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lattice is given by a product of complete elliptic integrals. It will be useful to consider the 
generating function D which has the explicit form 



:= ^^AVoJo^ = ' ^ (183) 

^V(i + 3y)(i-y)^ AGM(y(r^W+yF,7(r+W^) 

with a complete elliptic integral that is complementary to that in (|64p for D{y). For 
\y\ < ^, we have the expansion 

oo 

D{y) = HeunG (9, 3; 1, 1, 1, 1; 9y^) = ^ a^y^fc (i84) 

A:=0 

in terms of the hexagonal lattice integers in (jl77p . 

The full story, for the f.c.c, s.c, b.c.c. and diamond lattices, is then provided by the 
four identities 




2 



X 



V a. (t) (-x(l + x)(l + 9x))^ 
^ ((l-3x)(l + 3x))2^+i ^^^^^ 



fc=0 



V ' ^ , (187) 

^ ((l + x)3(l + 9x))^+- 



^ Q w X -I- I I 1 -I- I I 2 



((l + x)(l + 9x))^ 

which are valid for sufficiently small x and were obtained by simplification of formulae 
in 0611171 SH]. 

At bottom, all 4 results (|185p to (|188p originate from the first paper [16] of Wilfrid 
Norman Bailey's adroit series on infinite integrals involving Bessel functions, in much the 
same way that our proof of ([33]) resulted from his second paper [l^. In [I6i Eq. 8.1], 
Bailey showed that Jq°° Jfj,{at)J,^{bt)Jp{ct) dt is given by a product of 2F1 hypergeometric 
functions. In |461l47l H8]. Joyce used this result to obtain diamond and cubic lattice Green 
functions, in three spatial dimensions, from the square of the Green function (|183p for the 
hexagonal (or "honeycomb" ) lattice in two spatial dimensions, for which we have given a 
HeunG form in ()184p equivalent to that given in [42, Eq. 4]. With x > 0, the explicit form 
of E in (jl38p provides an analytic continuation for 

D ( = ^ ^ i^^) ■ (189) 

The singular value ki^, obtained from the diamond and f.c.c. lattice sums in (jl82p . 
also appears in the s.c. and b.c.c. lattice sums 

00 /2fc\ o /F 

E7Z|4«^ = ^^15^5/3 (190) 

fe=0 ^ ' 

^^24^(^/1+1)8^ " ^^'^ ^^^^^ 
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obtained by setting x = X15 := 3 — |\/5 k. 0.018576 in (jl86p and (jl87p . Setting x = X15 
in ()189p . we obtain the singular value ^5/3 = su\{a{2>y/x\^)) from definition of the angle 
a{x) in (jl39p . This is confirmed by the formula for tan(2a(x)) in p42p . We remark 
that (|190p follows from the functional relationship given in [501 Eq. 3.23] and that ()19ip 
follows from [24t Table 5.2a (A^ = 15)] and the Clausen product in [Ml Th. 5.7(a)(i), 
p. 180]. We now show how to obtain other singular values from lattice Green functions, 
using a modular identity. 

5.11 Cubic modular equations 

The cubic modular equation [24[ Th. 4.1, p. 110] 

OMHi) + O2iq)02iq^) = OMOsiq^) (192) 

relates instances of the Jacobi functions 62, in (jl54p . and ^3, in (jl03p . with nomas q and 
q^, to corresponding instances of 

00 

Vk^93{q) = eM--= E (-i)V (193) 

71= — 00 

with k' := Vl — fc^ and nome q associated with k. If we associate with the complete 
elhptic argument /, then (|192p gives the identity [211 Eqn. 4.2.6] \/kn/ + Vkl = 1. The 
results of the previous subsection follow from the notable circumstance that (|168p gives 

yjk'_^{z)k'_{z) + ^k+{z)k^{z) -1 = (194) 

which may be proven symbolically, by setting z = 1 — and denoting the left-hand side 
of (|194p by y{t). Then y{t) is analytic on the closed unit disk, \t\ < 1, and Maple computes 
an algebraic equation of the form y{t) P{y{t),t) = with P(0,0) non-zero. This proves 
that y{t) vanishes in some neighbourhood of t = and hence for |t| < 1. 

The resulting modular identities for the Green functions of cubic lattices are most con- 
veniently obtained from the Green function W4, for the diamond lattice, with a parametric 
solution 

W,{z,) - e,{q)9,{q ) , ^, (195) 

corresponding to the series solution ()188p at x = (1 — 2ry)/(3 + 6r]). Then the Green 
functions for the f.c.c, s.c. and b.c.c. lattices are given by the functional relationships 



^(1 - r/2) Wi{z,) = ^{l - r,^) W2{Z2) = J^—^Wsizs) = W^iz,) (196) 
3 6rj \ 1 + r] 



with the arguments 



^1 = 1-4??^, Z2 = \lzi{l-W zz = \—^\fzi, (197) 



r? V ' 2 + 2?? 

obtained from (jl85p to (|187p . respectively, by the substitution x = (1 — 2ry)/(3 + 6ry) that 
gave (|195p . We note the alternative b.c.c. parameterization [24^ Th. 5.7(a)(i), p. 180], 



W3(Z3) = 3F2 



111 

2' 2' 2 
1,1 
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The equivalence of the forms for 2:3 in ()197p and ()198p results from the modular identity 



(1 - 277)2 1 - 4r/2 



(2 + 27?)^ 



4/2(1-/2) = (2// 



/\2 



(199) 



with r] defined in (jl95p and / in (jl98p . To prove ()199p . we combined the cubic modular 
identity (fT92D with the Joubert-Cayley result [24, (4.6.14)] S{?>el{q^) /el{q)) = 0, where 
S is defined in (|156p . Alternatively, we may again apply the "modular machine" |251 §3]H 

Noting that the lattices sums in (jl8ip and (|182p yield the singular values /cs and fcis, 
from rational summands, we wondered if any other singular value might be obtained from 
a lattice Green function in such a neat manner. We know from Watson's classic work |71] 
that W3(l) yields the singular value k\ = l/\/2 while ^2(1) yields the singular value 
/cg = (V3 - \/2)(2 - \/3), as noted in [23]. Moreover, (fT98]) is equivalent to W:i{2kNk'j^) = 
AK'^/tt^, for = exp(— 7r\/iV)- Hence ^3(^3) yields the singular values ki, k^, kj for the 
rational arguments 23 = 1, i, |, respectively. 

Prompted by the sum over s.c. lattice integers in (|190p . we sought further examples, 
in which a sum over rational numbers might lead to a singular value /cat. We found 5 new 
results for W2, which appear to exhaust the cases with rational summands. These occur 
with A^/3 = 7, 11, 19, 31, 59, for which we obtained 
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\/3-l) y/bOiKmK, 



1 77 -f^ 59/3 
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(3757-13^/3) K57i^i9/3 = -^G{b7) 
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— 39^3 - 7^/93 K93i^3i/3 = TTT G(93) 



31 

23^/3 



(200) 



— (3733 - 57TT) K33i^ii/3 = ^ G(33) (201) 



(202) 
(203) 



G(177) (204) 



with reductions to T values given by 



G{N) 



2y27r n 

^ n=l 



\(-AN\n) 



(205) 



where T{n/ (4A^)) contributes to the product if n is coprime to AN and then occurs with an 
exponent ib|, according as the sign of the Legendre-Jacobi-Kronecker symbol (— 4A^|n). 
We remark that p = 7,11,19,31,59 are the only primes for which = 3p is a disjoint 
discriminant of type one, as considered in [24i Eq. 9.2.8, p. 293]. 



^Striking cubic modular equations [q^ oy N 9A^), originating with Ramanujan, are explored in [24l 



§4.7]. In particular there are attractive cubic recursions for the cu6ic multiplier M = \/(l 
9Ul^)/dU'l)^ as occurs in (fTQej) . 
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We may obtain other singular values by choosing the argument of W4 in (jl95p to be 
an appropriate algebraic number. For example, the sums 



V = 4^39^^13/3 (206) 



k 



g(i7fT^)^WfTI7fP = ^^105i^35/3 (207) 

over the diamond lattice integers bk in (|177p have relatively simple surds in their sum- 
mands, obtained by setting q = exp(— tty^IS/S) and q = exp(— 7ry^35/3) in the parametric 
solution (jl95p . We may also obtain quartic values of Z4 and evaluations like (j207p at the 
even singular values with A^/6 = 3,5,7,13,17. We remark that p = 5,7,13,17 are 
the only primes for which N/2 = 3p is a disjoint discriminant of type two, as considered 
in PI Eq. 9.2.9, p. 293]. 

5.12 Integral sum rules 

Unfortunately, the discrete sum rule (11811) . at the singular value ks, does not prove con- 
jecture (jlOip but instead converts it to the conjectural integral sum rule 

thit) {lo{t)Ko{2t) - ^ Io{2t)Ko{t)j K^{t) dt =1 (208) 

where the term containing Io{2t) is now proven to yield the singular value k^, whose 
appearance in the term containing KQ{2t) was conjectured in (jlOip . 
It looks to be an even tougher proposition to prove the sum rule 

l^^tloit) (^Ko{t) - ^/o(t)) Kl{t)dt =1 (209) 

for which we have now obtained two representations in terms of integrals of complete 
elliptic integrals, namely 

/^^£kLfarctanhf,/IH|^V^V^ = « W 



farcsinhP)-^V- = « (211) 



V4 + x2 V ^2/ ^ 

with D{y) given in E{x) given in (|138p and the 3-loop sunrise diagram 55^1 appearing 
via the more demanding logarithmic terms in (jl23p and (jl40p . 

As a companion to the proven evaluation ([50|) and the discrete sum rule ([52]) . we 
present 

farcsinh(2tan0) - ^) dO (212) 

io Vl + Ssin^^ V ' ' 

as our penultimate conjecture, also checked to 1200 decimal places. If the mathematics 
of sum rule (j212p might be elucidated at the singular value A;3, then there might be some 
hope for a proof of the quantum field theory result (|209p , at the singular value /cis . 
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5.13 The even moment C5 







We were unable to derive single integrals of elliptic integrals for the even moments c^^2k- 
A double-integral representation is readily available by setting w = in the Aufbau 



SQ{a,b,c,d,e,w) = j S4{a,b,c,v)S3{d,e,v + w) dv 

83(0, b, u)52(c, u + v)S3{d, e,v + w) du dv (213) 



—00 

00 fOO 



00 J —00 



to obtain _ _ 

2 J -00 J -00 + (n + vy 

with 53 given by ()30p . Setting the 5 parameters to unity and making the transformations 
u = 2 tan 6 and v = 2 tan (j), we obtain 

„ = = r r '^'""''>'^'""t> do . (215) 

2 J_^/2 J-7r/2 V COS^ 6* COS^ (/> + 4 siu^ (6' + (/)) 

The higher even moments are obtained by suitable differentiations of (j214p with respect 
to c, before setting c = 1. 



6 Six Bessel functions 

We are now equipped to write the odd moments te.i) -se,! and cg,! as single or dou- 
ble integrals over complete elliptic integrals, with integrands that are computable with 
great efficiency, using the exponentially fast process of the arithmetic-geometric mean, 
discovered by Lagrange, around 1784, and independently by Gauss, at the age of 14, in 
1791 [3l[24]. 

6.1 The odd moments t6,2fc+i 

We begin by folding D4, in (f62]) . with T4, in (I147p . to determine 

Te{u,a,b,c,d,v) := / tjQ{ut)KQ{at)Ko{bt)Ko{ct)KQ{dt)Jo{vt)dt 

Jo 

/•oo 

= / 2w D4{a,b,c,w)T4{u,d,w,v) dw 

Ja+b+c 

where we group the 3 internal lines with masses a, b and c to have a total momentum with 
norm w'^. Setting a = b = c = d=l, u = v = i and w = 1/y, we obtain one instance of 
D(y), from its definition D{y) := 21)4(1, 1, 1, l/y)/y, and another, somewhat surprisingly, 
from the novel result D{y) = 6T4(i, 1, l/y,i)/y in (jl49p . Hence we obtain 

i6,i:= rtll{t)Kl{^^t = \ r^dy (216) 

Jo -J ^0 

with an initial factor of | arising from (jl37p . in the special case (jl35p . This is a rather 
efficient representation of te.i) which delivers 1200 decimal places in two minutes, using 
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Pari-GP. Similarly, we may derive a one-dimensional integral for higher odd moments 
by differentiations of D4{a,b,c,w) with respect to one of its masses, before going to the 
equal-mass point. More conveniently, one may use the summation 



n=0 ^ ^ 



2n\ C49fc- 



'2n+l 



(2"n!)2 



that follows from ([^^j) . By this means, we found that H^e.s — reproduces the value 

of t6,5 to 1200 decimal places. 

6.2 The odd moments S6,2/c+i 

We were unable to derive a one-dimensional integral over agm functions for sg^i, though 
we shall conjecture such an integral, in the next subsection. Here, the best that we can 
do comes from using the folding 

POO 

SQ{a,b,c,d,e,w) = / 2u D4{a,b,c,u)S4{d,e,u,w) du 

J a+b+c 

which leads to a choice of integrals for 

POO poo 

S4{d,e,u,w) = / 2vD4{d,e,u,v)S2{v,w)dv= / 2vDs{d,e,v)S3{v,u,w)dv 

J d+e+u J d+e 

with the first form involving an agm and the second an arctanh procedure. The latter is 
more convenient, since it yields the rectangular double integral 



^-'"f'"''^™ d= (217) 

Jo Jo Y^(l - 42;2)A+A_ 

with A± := (y it z)"^ — y'^z'^, obtained by transforming to y = 1/u and z = 1/v, and D{y) 
given in (j64p . We remark that is positive, within the rectangle of integration, and that 
when A- is negative an analytic continuation of arctanh to arctan keeps the integrand 
real and positive. Higher odd moments may be obtained by differentiations with respect 
to w, before setting w"^ = —1. More conveniently, we may use the summation 



oo 



n=0 ^ '' 

since cs^i and 05^3 were computed by Broadhurst to 200,000 decimal places (see footnote 
3) and higher moments are (conjecturally) determined by them, using (jl20p and the 
appropriate recursion from ([TO]) . By this means, we found that ||s6,3 — ^^6,1 reproduces 
the value of S6,5 to 1200 decimal places. 

6.3 Sum rules 

Our final conjecture is that there is a infinite tower of sum rules relating moments in 
which powers of Kq are replaced by corresponding powers of ttIq. We were alerted to this 
possibility by the sum rule 

poo 

/ {iT^ll{t)-Kl{t))Kl{t)dt = Q (219) 
Jo 
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proven in (j9ip . For each pair of integers (n, k) with n >2k > 2 we conjecture that 

^"^^^ / n \ r°° 7fiil 

Z2n,n-2k:= E^-^rio / [^^0(0]""'™ [i^o dt (220) 

m=0 ^ 



with the vanishing of Z^^q proven in (j219p . in the case with n = 2 and k = 1. 
With 6 Bessel functions in play, the sum rule 

Ze,i := TT / tIo{t) (7r2/2(t) _ ^K^^t)) K^{t) dt (221) 
Jo 

relates a pair of odd moments, giving the conjectural evaluation 

56,1 = ^ rtll{t)K^{t)dt. (222) 

-J JO 

At first sight, this seems to be hard to check, at high precision, because it involves the 
slowly convergent moment 



t/o^(t)Ko^(t)dt = ^a,|gtl 



k=0 



with an integrand of order l/t"^, at large t, and a summand of order at large k, 

coming from the hexagonal lattice integers in (j24|) . However, we were able to exploit 
the integral representation 

1 



to obtain 

tr^it)K^{t) dt= r D{y)D{y) dy (223) 



as an integral over a pair of agm functions, with D[y) given by ()64p and D{y) by (jl83p . 
Then conjecture (j22ip is equivalent to the evaluation 

56,1 := r th{t)Kl{t)dt [' D{y)D{y)dy (224) 

Jo "J Jo 

which was confirmed to 1200 decimal places, by setting A; = in (|218p to compute sg^i 
and by using Pari-GP to evaluate the integral over y. We find it remarkable that the 
complicated double integral 4-loop sunrise diagram in (j217p seems to be reducible to the 
attractive single integral in (j224p . by removing a factor ^vr^. 

With 8 Bessel functions in play, conjecture (|220p gives a pair of sum rules. Prom the 
vanishing of Z^^ and .^8,2, we obtain 

C8,o = r ll{t)[QKl{t)-^^ll{t))K^{t)dt (225) 
Jo 

C8,2 = r t'll{t){QKl{t)-^''ll{t))Kl{t)dt (226) 

Jo 
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again noting the slow convergence of these integrals and of the equivalent sums 

C8,o = YAi^) ( f^,^ )6c6,2fc -vr26fcC4,2fc ) (227) 



k=0 ^ ^ 

oo 



?[ii] / vr \2 f/2k\ 2r \ / N 

C8,2 = 2^ ) [[ k )^'^6,2fc+2 - vr 6fcC4,2fc+2 1 (228) 

k=0 ■ \\ / / 

which involve the diamond lattice integers bk in (j56p or p77p . 

Conjecture (|220p gives novel evaluations of C4n,2k and S4n+2,2A:+i5 for n > /c > 0. We 
remark, however, that p27p and ()228p do not exhaust the integer relations for moments 
with 8 Bessel functions. We also found that the ratio 

coincides with ^vr^, to 80 decimal places. 

6.4 The odd moments C6,2fc+i 

Grouping the 6 internal lines of the 5-loop vacuum diagram 

/"OO fOO 

VQ{a,b,c,d,e, f) = / / 4uv D4{a,b,c,u)V2{u,v)D4{v,d,e, f) dv du 

Ja+b+c Jd+e+f 

into two sets of 3 lines, we obtain 

fi„^^ fimiog^ (229) 

Jo Jo z - y 

after setting the masses to unity and transforming to y = 1/u and z = 1/v. We also note 
that Broadhurst conjectured that the value of cg^s is ||c6,3 — ^C6,i + This was later 
checked to 500 decimal places, using data in [9]. 

6.5 The even moment cg o 

Finally, we obtain 3 complete elliptic integrals in the integrand of 

/oo poo 
/ S3{a,b,u)S3{c,d,u + v)S3{e,f,v)dudv (230) 
-oo J ~oo 

using the Aufbau (|29p . This then delivers 

,^/2 .^/2 K(sin 6) K(sin </.) K f E2M= 



^ J-7r/2 J-7r/2 VCOS^ 9 COS"' (/> + Sm^ [9 + (j)) 

by the same transformations as for ()215p . 
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7 Computational notes 



This paper contains several proofs of identities that we first conjectured on the basis 
of numerical investigation, hugely facilitated by access to Sloane's wonderful sequence 
finder. For the many one-dimensional integrals that we have noted, we were greatly 
aided by the efficiency of the agm and intnum procedures of Pari-GP, for evaluations of 
integrands and integrals at precisions up to 1200 decimal places. Results were then fed to 
the implementation of the PSLQ algorithm in Pari-GP^s lindep procedure, with which 
we performed many unsuccessful searches for integer relations, as well as obtaining the 
positive results reported in the paper. By way of example, we remark that the integral 



was evaluated to high precision in order to search for relations between cs^i, 05^3, £ and 
products of powers of vr, C and 1/C, with coefficients that might be Q-linear combinations 
of 1, \/3, \/5 and \/TE. No such relation was found. 

Maple was especially useful for its HeunG, MeijerG and Gosper procedures and also 
for quick PSLQ searches with few terms, at relatively low precision. 

But for two-dimensional numerical quadratures we found neither Pari-GP nor Maple 
to be remotely adequate for our demanding investigations. For these, we came to rely 
on Bailey's multiprecision code for two-dimensional integrals 111', '14], which confirmed, 
to more than 100 decimal places, the correctness of derivations of (j215p . (j23ip and other 
identities. Here we offer a brief description of this scheme. 

7.1 Multi-dimensional quadrature 

Bailey's 1-D and 2-D schemes, as well as the one-dimensional intnum procedure in Pari- 
GP, employ the tanh-sinh quadrature algorithm, which was originally discovered by Taka- 
hasi and Mori [67j . It is rooted in the Euler-Maclaurin summation formula [41, p. 285], 
which implies that for certain bell-shaped integrands /(x) on [0, 1] where the function and 
all higher derivatives rapidly approach zero at the endpoints, approximating the integral 
of f{x) by a simple step-function summation is remarkably accurate. This observation 
is combined with the transformation x = tanh((7r/2) sinht), which converts most "rea- 
sonable" integrand functions on ( — 1,1) (including many functions with singularities or 
infinite derivatives at one or both endpoints) into bell-shaped functions with the desired 
property. 

In particular, we can write, for an interval length h > 0, 



where g{t) = tanh((7r/2) sinh t), Xj = g{hj), Wj = g'{hj), and N is chosen large enough 
that \wjf{xj)\ < e for \j\ > N. Here e = 10"^, where p is the numeric precision level in 
digits. Note that the resulting quadrature formula (I232p has the form similar to Gaussian 
quadrature, namely a simple summation with abscissas xj and weights wj, both sets of 
which can be pre-computed since they are independent of the integrand function. For 
many integrand functions, once h is sufficiently small, reducing h by half yields twice as 
many correct digits in the result (although all computations must be performed to at 





(232) 
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least the level of precision desired for the final result, and perhaps double this level if the 
function is not well-behaved at endpoints). Additional details of efficient implementations 
are given in [SJ [T3j . 

One of the numerous applications of one-dimensional tanh-sinh quadrature in this 
study was the verification of our final conjecture given in (j22U|) . This was done using 
Bailey's implementation of the one-dimensional tanh-sinh algorithm, together with the 
ARPREC extreme-precision software package |13j . Evaluation of the Bessel function Io{t) 
was performed using a hybrid scheme where formula 9.6.12 of [T| was used for modest- 
sized values, and formula 9.7.1 for large values. Evaluation of i^o(^) was performed using 
9.6.13 of [U for modest-sized values, and 9.7.2 for large values. Note however that formula 
9.6.13 for Ko^t) must be implemented using a working precision that is roughly twice the 
level desired for the final result, due to the sensitive subtraction operation in this formula. 
Also note that when m = in (I220p . this combination of formulas is not satisfactory, 
because for large t the function /o(*) is exceedingly large, and Ko{t) is exceedingly small, 
and even though the product is of modest size, overflows and underflows are possible in 
intermediate function evaluations, even when using high-precision software that has an 
enormously extended dynamic range. For such cases (m = and t large), we employed 
formula 9.7.5 of [l], which gives an asymptotic series for the product Io{t)KQ{t). We had 
differently addressed this issue in the special case of (j222p . 

Armed with an efficient implementation of these schemes, we were able to verify (I220p 
for all (n, k) pairs, where 1 < fc < [n/2j and 4 < n < 12 (there are 43 such pairs), in each 
case to over 340-digit accuracy. In addition, we verified ()220p for a variety of larger (n, k) 
pairs, including (15,7), (20,10), (25,11), (30,12), (37,13), and (41, 14), again to over 
340-digit accuracy in each case. 

The tanh-sinh quadrature algorithm can also be performed in two or more dimensions 
as an iterated version of the one-dimensional scheme. Such computations are many times 
more expensive than in one dimension. For example, if roughly 1,000 function evaluations 
are required in one dimension to achieve a desired precision level, then at least 1,000,000 
function evaluations are typically required in two dimensions, and 1,000,000,000 in three 
dimensions. Additionally, the behaviour of multi-dimensional tanh-sinh quadrature on 
integrand functions with singularities or infinite derivatives on the boundaries of the region 
is not as predictable or well-understood as in one dimension. 

Nonetheless, we were able to use 2-D tanh-sinh quadrature to successfully evaluate 
a number of the double integrals mentioned in this paper, after making some minor 
transformations. As one example, consider the integral mentioned above for c^^, namely 



Note that this function (see Figure 1) has singularities on all four sides of the region of 
integration, with particularly troublesome singularities at (7r/2, — 7r/2) and (— 7r/2, 7r/2). 
However, after making the substitutions s <^ tt/2 — s, t <— tt/2 — t and r ^ s/t, and 
taking advantage of the symmetry evident from Figure 1, we obtain 




(233) 




Jo Y^sin^(rt) sin^ t -|- 4sin^(t(l -|- r)) 
1-77/2 H tK(cos(rt))K(cost)drdt 




(234) 
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Figure 1: Plot of C5 integrand function in fl233p 



which is significantly better behaved (although these integrands still have singularities 
on two of the four sides of the region). As a result, we were able to compute cs^o with 
this formula to 120-digit accuracy, using 240-digit working precision. This run required 
a parallel computation (using the MPI parallel programming library) of 43.2 minutes 
on 512 CPUs (1024 cores) of the "Franklin" system at the National Energy Research 
Scientific Computing Center at the Lawrence Berkeley National Laboratory. The final 
result matched the value that we had previously calculated using ([5]) (see [9j) to 120-digit 
accuracy. 

This same strategy was successful for several other 2-D integrals. For example, we 
computed C6,o to 116-digit accuracy, which again matched the value we had previously 
computed, in 64.8 minutes on 1024 cores of the Franklin system. In the case of C6,i, the 
transformation described above for cs^o converted the integrand function of (I229p into a 
completely well-behaved function, without any singularities. As a result, we were able to 
compute C6,i to 120-digit accuracy using only an Apple Intel-based workstation with four 
computational cores, in 28 minutes. As before, the result matched the earlier calculation. 

As already noted, complex numbers are avoided in integral (|217p by writing it as 




(235) 



where 




for X > 
for X = 
for X < 
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yields positive real numbers within the rectangle of integration. We were able to confirm 
that the double integral (j235p yields the first 120 of the 1200 decimal places obtained, far 
more easily, from the single integral (|224p . 

8 Conclusions 

Despite some notable progress in discovering and proving new results, we are left with 8 
outstanding conjectures]! Of these, 5 have their first instances in Equations ([Mj) to ([97|) 
and Equation (flOT]) of Section 5.1, with 3 outliers, in Equations (fT20]) . (|2T2]) and (12201) . 

Conjecture (f94j) lies deep in 4-dimensional quantum field theory, but it is reasonable 
to suppose that it might be derivable from the two-dimensional conjectures (f95|) . ([96]) 
and (I97p . together with their proven sign-changed variants in (IDSh . (I99p and (llOOp . 

The conjectural integer relations ()120p and ()125|) may be provable by adding rational 
data such as (jl2ip to a set of recursions richer than those considered in |26] [60] . 

The real challenge is set by the remarkable sum rule ()209p . with a dispersive presen- 
tation (|210p . a non-dispersive presentation ()21ip and a kindergarten analogue (|212p . 

Acknowledgements. We are most grateful to Andrei Davydychev, for elucidation 
of dilogarithms obtained in quantum field theory, to Stefano Laporta, for the astute 
observations in 4-dimensional quantum field theory from which our work arose, to Bas 
Tausk, for advice on Mellin-Barnes transforms, and to John Zucker, for several pertinent 
observations regarding integrals of elliptic integrals. 



The conjecture in pOip was later proven by contour integration — see [27] . 
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